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INTRODUCTION 

This report covers our research activities from October 25, 1983 to 
October 24, 1984 under the auspices of NASA Grant No. NAG 3-477; the 
report is actually a continuation of our previous report^, NASA CR- 
174661, for NASA Contract NAS 3-23542. In the previous work, soot 
formation in toluene-, benzene-, and acetylene-oxygen-argon mixtures was 
investigated. We have continued the study of the influence of oxygen on 
soot formation from hydrocarbons; this report presents the results 
obtained in oxidation of allene, 1,3-butadiene, vinylacetylene, and 
chlorobenzene . 
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We have suggested ’ that the soot formation mechanism is probably 
the same for both pyrolysis and oxidation of hydrocarbons. Therefore, 
in continuation of our efforts to investigate the effects of fuel 
structure on the chemical kinetic mechanisms responsible for soot 
formation, we have concentrated on the pyrolysis studies. We report 
here our new results obtained in shock-tube pyrolysis of ethylene, 
vinylacetylene, 1-butene and chlorobenzene. We have initiated a study 
of soot formation in mixtures of hydrocarbons. The report presents the 
results obtained in pyrolysis of acetylene-hydrogen, benzene-acetylene, 
benzene-butadiene, and chlorobenzene-acetylene argon-diluted mixtures. 

The report is organized in three main pa ’ts: experimental results, 

discussion of the experimental results, and computer modeling of soot 

formation in shock- tube pyrolysis of acetylene. The latter is the first 

attempt of a detailed chemical kinetic modeling of soot formation. The 

3 

modeling work was performed in collaboration with Professor William C. 
Gardiner, Jr. of the University of Texas and Dr. Stephen E. Stein of the 
National Bureau of Standards. 
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APPARATUS AND PROCEDURES 

The experimental apparatus and procedures used in this study were 
similar to those described in our previous report*. Briefly, the 
experiments were conducted behind reflected shock waves in a 7.62 cm 
i.d. shock tube. The test gas mixtures were prepared manometrically. 
The stated purities of the gases were; argon-99 . 998%, acetylene-99. 6%, 
oxygen-99 .3%, allene-93%, 1 ,3-butadiene-99%, ethylene-99 .5%, 1-butene- 
99.9% and hydrogen-99 . 99%. These gases were used without further 
purificatior The vinylacetylene used was of technical grade with up to 
2% impurities. It was supplied by Wiley Organics as a 50% solution by 
weight in xylene, to which 0.1% by weight of P-tertiary-butylcatechol 
had been added. From this solution vinylacetylene vapor was distilled 
at 0°C and was found by gas chromatographic analysis to contain less 
than 0.5% xylene. The toluene (Reagent, Baker), benzene (Spectrana- 
lyzed, Fisher) and chlorobenzene (Reagent, Baker) were purified by 
repeated freezing and evacuation. The shock tube was cleaned after 
every run. 

The soot conversion was det* rmined by measuring the attenuation of 
a He-Ne 632.3 nm laser beam. As in our previous studies, the soot 
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yields were calculated according to Graham's model but leaving out the 

2 2 

quantity E(m) = - Im[(m -l)/(m +2)], where m is the complex refractive 
index of soot particles. This arbitrary form of reporting the results 
was chosen to emphasize the ambiguity in the value of m and in the 
laser-extinction model itself^ ' . An additional source of uncertainty 
is possible absorption of the laser beam by molecular species rather 

than soot particles; light at 632.8 nm is already absorbed by compounds 

7 8 9 

with as few as six aromatic rings . However, Rawlins and co-workers ’ 
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concluded i>om their recent optical studies thct light extinction at 
632,8 nm is primarily caused by soot particles. 

TEST RESULTS 

Table I summarizes the experimental conditions of the mixtures 
covered in this report. Mixtures A to K were used to study and compare 
soot formation from aliphatic hydrocarbons: acetylene, ethylene, 

vinylacetylene and 1-butene. T.le remaining mixtures (L to W) were 
utilized to study soot formation from individual aromatic hydrocarbons 
(benzene and chlorobenzene) and their mixtures with aliphatic 

hydrocarbons (acetylene, allene, vinylacetylene and 1 ,3-butadiene) . The 
measured soot yield, induction times for soot appearance and rates of 
soot formation are reported in Figs. 1 through 51. 
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TABLE I. EXPERIMENTAL CONDITIONS. 


Mixture 

Composition 
(% vol. in argon) 
Fuel O 2 

T 5 

(K) 

P 5 

(bar) 

C 5 x 10 6 
(mol/cm 3 ) 

(C) x 10“ 17 
(carbon- 
atoms/ cm J ) 

A* 

1.09, C 2 H 2 

1665-3040 

2.14-3.87 

15.31-15.69 

2.01-2.06 

B 

2.72 , C 2 H 2 

1674-2413 

2.11-3.05 

15.16-15.42 

4.97-5.05 

C* 

4.65, C 2 H 2 

1687-3122 

1.25-2.33 

8.70- 9.22 

4.87-5.16 

D 

4.65, C 2 H 4 

2024-2874 

1.50-2.11 

8.68- 8.98 

4.86-5.03 

E* 

0.726,C 3 H 4 

*'26-2294 

1.93-2.91 

15.06-15.39 

1.98-2.02 

F 

0.726,C 3 H 4 0.726 

1497-2041 

1.93-2.56 

15.10-15.72 

1.98-2.06 

G* 

0.54, C 4 H 6 

1664-2415 

2.10-3.01 

14.94-15.32 

94-1.99 

H 

0.54, C 4 H 6 0.54 

1560-2082 

2.02-2.68 

15.09-15.65 

1.96-2.04 

I 

0.54, C 4 H 4 

1670-2424 

2.18-3.09 

15.31-15.73 

1.99-2.05 

J 

0.54, C 4 H 4 0.54 

1599-2020 

2.08-2.60 

15.40-15.71 

2.00-2.04 

K 

0.54, C 4 H 8 

1630-2299 

2.08-2.96 

15.34-15.66 

2.00-2.04 

L + 

0.311, C ? Hg 

15252349 

1.98-3. CO 

14.52-15.74 

1.90-2.06 

M + 

0.311, CgHg 

1561-2272 

1.93-2-87 

14.90-15.59 

1.67-1.75 

N + 

0.311, C 6 H 6 0.311 

1515-2242 

1.98-2.83 

15.16- 15.68 

1.70-1.76 

P 

0.311, C 6 H 5 Cl 

1446-2201 

1.89-2.81 

15.1i-15. 86 

1.70-1.78 

Q 

0.311, C 6 H 5 CI 0.311 

1433-2184 

1.81-2.77 

14.74-15.74 

1.66-1.77 

R 

4.65, C 2 H 2 
+ 4.65, H 2 

1744-2427 

1.30-1.80 

8.84- 9.08 

4.. 5-5. 09 

S 

0.311, C 6 H 5 C1 
+ 1.09, C 2 H 2 

1475-2237 

1.92-2.82 

15.02-15.75 

3.66-3.84 

T 

0.311, CgHg 
+ 0.54, C 2 H 2 

1571-2300 

2.03-2.91 

14.53-15.58 

2.58-2.76 
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Table I. 

(continues) 



Mixture 

Composition 
(% vol. in argon) 
Fuel 0 2 

T 5 

(K) 

P 5 

(bar) 

C 5 x 10 6 
3 

(mol/cm ) 

[ C ) _ X io -1 
(carbon^ 
atoms/cin ) 

U 

0.311, C 6 H 6 
+ 1.09, C 2 H 2 

1506-2279 

1.94-2.86 

15.10-15.52 

3.68-3.78 

V 

0.311, C 6 H 6 
+ 2.72, C 2 H 2 

1593-2427 

2.04-3.06 

14.97-15.55 

6.58-6.84 

W 

0.311, C 6 H 6 
+ 0-54, C 4 H 6 

1537-2433 

1.98-3.07 

15.14-15.51 

3.67-3.76 


* These experiments were performed with the partial support of the Office of Fossil 
Energy, U.S. Department of Energj , under the auspices of Grant Number DE-FG 22-80 
PC30247. 


These experiments were reported in NASA CR-174661 and are reproduced here for 
purposes of comparison. 
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Figure 2. Comparison of soot yield at 1.0 ms for the mixtures 
2.72% acetylene-argon and 1.09% acetylene-argon. 
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Figure 4. Soot formation rates for the mixture 2.72% 
acetylene-argon. 
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Figure 5. Soot yields vs. temperature at different reaction 
times for the mixture 4.65% ethylene-argon. 
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Figure 6. Comparison of soot yield at 1.0 ms for the mixtures 
4.65% acetylene-argon and 4.65% ethylene-argon. 
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Figure 7. Soot induction times for the mixtures 4.65% 
ethylene-argon and 4.65% acetylene-argon. 


7.00 



13 



Figure 8. Soot formation rates for the mixtures 4.65% 
ethylene- argon and 4.65% acetylene-argon. 
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Figure 9. Soot yields vs. temperature at different reaction 
times for the mixture 0.726% allene-0. 726% oxygen- 
argon . 
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Figure 10. Comparison of oot yield at 1.0 ms for the mixtures 
0.726% allene-argon and 0.726% allene-0. 726% 
oxygen-argon. 
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Figure 11. Soot induction times for the mixtures of 0.726% 

allene-argon and 0.726% allene-0. 726% oxygen-argon. 
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Figure 12. Soot formation rates for the mixtures 0.726% 
allene-argon and 0.726% allene-0 . 726% oxygen- 
argon. 
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Figure 13. Soot yields vs. temperature at different reaction 
times for the mixture 0.54% 1 ,3-butadiene-0 .54% 
oxygen-argon. 
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Figure 14. Comparison of soot yield at 1.0 ms for the mixtures 
0.54% 1,3-butadiene-argon and 0.54% 1 ,3-butadiene- 
0.54% oxygen-argon. 
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Figure 15. Soot induction times for the mixtures 0.54% 1,3- 
butadiene-argoa and 0.54% 1 ,3-butadiene-0.54% 


oxygen-argon. 
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Figure 16. Soot formation rates for the mixtures 0.54% 1,3- 
butadiene-argon and 0.54% 1 ,3-butadiene-0 . 54% 
oxygen-argon. 
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Figure 17. Soot yields vs. temperature at different reaction 
times for the mixture 0.54% vinylacetylene-argon. 



ld m xE cm) 


23 


D 

c n 



TEMPERATURE (K) 


Figure IS. Soot yields vs. temperature at differeut reaction 
times for the mixture 0.54% vinylacetylene-0 .54% 
oxygen-argon. 
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Figure 19. Comparison of soot yield at 1.0 ms for the mixtures 
0.54% vinylacetylene-argon and 0.54% vinyla cetyl ene- 
0.54% oxygen-argon. 
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Figure 20. Comparison of soot yield at reaction times of 0.5, 
1.0, 1.5 and 2.0 ms for the mixtures 0.54% vinyl- 
acetylene-argon and 0.54% vinylacetylene-0.54% 
oxygen- argon. 
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Figure 21. Soot induction times for the mixtures 0.54% 

vinylacetylene -argon and 0.54% vinylacetylene-0 .54% 
oxygen-argon. 
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Figure 2J. Soot yields vs. temperature at different reaction 
times for the mixture 0.54% 1-butene-argon. 
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Figure 25. Soot induction times for the mixtures 0.54% vinyl 
acetylene-argon, 0.54% 1-butene-argon and 0.54% 
1,3-butadiene-argon. 
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Figure 26. Soot formation rates for the mixtures 0.54% 

vinylacetylene-argon, 0.54% 1-butene-argon and 
0.54% 1,3-butadiene-argon. 
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Figure 28. Comparison of soot yield at a reaction time of 1.0 
ms for the mixtures 0.311% chlorobenzene-argon 
and 0.311% benzene-argon. 
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Figure 29. Soot induction times for the mixtures 0.311% 
benzene-argon and 0.311% chlorobenzene-argon. 
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Figure 30. Soot formation rates for the mixtures 0.311% 
benzene-argon and 0.311% chlorobenzene-argon. 
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Figure 31. Soot yields vs. temperature at different reaction 
times for the mixture 0.311% chlorobenzene-0. 311% 


oxygen- argon. 
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Figure 32. Comparison of soot yield at a reaction time of 1.0 
ms for the mixtures 0.311% chlorobenzene-argon and 
0.311% chlorobenzene-0. 311% oxygen-argon. 
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Figure 33. Comparison of soot yield at a reaction time of 

1.0 ms for the mixtures 0.311% chlorobenzene-0. 311% 
oxygen-argon and 0.311% benzene-0 .31 1% oxygen-argon. 
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Figure 34. Soot induction times for the mixtures 0.311% 

chlorobenzene-argon and 0.311% chlorobenzene-0 . 311% 
oxygen-argon. 
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Figure 35. Soot formation rates for the mixtures 0.311% 
chlorobenzene-argon and 0.311% chlorobenzene- 
0.311% oxygen- a rgor.. 
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Figure 36. Soot yields vs. temperature at different reaction 
times for the mixture 4.65% ?.cetylene-4.65% 
hydrogen-argon . 
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Figure 37. Comparison of soot yield at a reaction time of 1.0 
ms for the mixtures 4.65% acetylene-argon and 4.65% 
acetylene-4 .65% hydrogen-argon. 


SOOT YIELD (X) xE (ITU 


A3 



TEMPERATURE (K) 


Figure 38. Soot yields vs. temperature at different reaction 
times for the mixture 0.311% chlorobenzene- 1 .09% 
acetylene-argon. 
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Figure 39. Comparison of soot yield at a reaction time of 
1.0 ms for the mixture*. 0.311% chlorobenzene- 
argon, 0.311% chlorobenzene- 1 .09% acetylene- 
argon and 1.09% acetylene-argon. 
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Figure 39a. Comparison of the amount of soot formed at a reaction 
time of 1.0 ms for the mixtures 0.311% chlorobenzene- 
argon, 0.311% chlorobenzene- 1 .09% acetylene-argon and 
1.09% acetylene-argon. 
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Figure 40. 


Soot induction times for the mixtures 0.311% 
chlorobenzene- 1 .09% acetylene-argon and 0.311% 
chlorobenzene-argon. 
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Figure 41. Soot formation rates for the mixtures 0.311% 

chlorobenzene- 1 .09% acetylene-argon and 0.311% 
chlorobenzene-argon. 
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Figure 42. Soot yields vs. temperature at different reaction 

times for the mixture 0.311% benzene-0. 54% acetylene- 


argon. 
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Figure 43. Soot yields vs. temperature at different reaction 
times for the mixture 0.311% benzene-1 .09% 
acetylene-argon . 
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Figure 44. Soot yields vs. temperature at different reaction 
times for the mixture 0.311% benzene-2, 72% 
acetylene-argon. 



51 


g T I M E = 1 . 0 ms 

»- RBS0RPTI ON RT 632.8 niti 

□ 0.311 ZC B H B +0. 54 V, C 2 H 2 +Rr 

CCD*w2. 71 nIO 17 carbon atoms/om 9 


g 0 0.311ZC B H B +2.72 Z C 2 H 2 +Rr 

in- CCD-^6. 71 kIO 17 carbon atoms/cm 3 



o 

o 


c\n 


A 0.311 ZC B H B + 1.09 XCgHg+fir 

CCZK3. 71 MlO 17 carbon atoms/cm 3 

+ 0.311 ZC B H„ + Rr 

CCJ-vl . 7 1 *10 17 carbon atoms/cm 3 



TEMPERRTURE (K) 


Figure 45. Comparison of soot yield at a reaction time of 1.0 
ms for the mixtures 0.311% benzene-0 . 54% acetylene- 
argon, 0.311% benzene-2. 72% acetylene-argon, 0.311% 
benzene-1 .09% acetylene-argon and 0.311% benzene-argon. 
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Figure 45a. Comparison of the amount of soot formed at a reaction 
time of 1.0 ms for the mixtures 0.311% benzene-0 .54% 
acetylene-argon, 0.311% benzene-2. 72% acetylene- 
argon, 0.311% benzene-1 . 09% acetylene argon and 0.311% 
benzene-argon. 
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Figure 46. Soot induction times for the mixtures 0.311% 

benzene-argon, 0 . 31 1% benzene-0 .54% acetylene- 
argon, 0.311% benzene-1 .09% acetylene-argon and 
0.311% benzene-2. 72% acetylene-argon. 




Logi 0 (Rate) , (sec -1 ) 



Figure 47. Soot formation rates for the mixtures 0.311% 

benzene -argon, 0.311% benzene-0 .54% acetylene- 
argon, 0.311% benzene-1 .09% acetylene-argon and 
0.311% benzene-2. 72% acetylene-argon. 


HBS0RPTION RT 632.8 nm 
0.311 X C s M.j + 0.54 ZC 4 M B + Rr 
CCD~3. 71 mIO 17 carbon atoms/cm 3 
P=1 .98 to 3.07 bar 

□ 0.5 ms 

O 1.0 ms /^\ 

A 1.5ms T \ 

+ 2.0 ms / \ 



1700. 1900. 2100. 

TEMPERRTURE (K) 


Figure 48. Soot yields vs. temperature at different reactions 
times for the mixture 0.311% benzene-0 . 54% 

1 , 3-butadiene-argon . 
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Figure 49. Comparison of soot yield at a reaction time of 1.0 
tns for the mixtures 0.311% benzene-0 .54% 
1,3-butadiene-argon, 0.311% benzene-argon, 0.311% 
benzene- 1 .09% acetylene-argon and 0.54% 1,3-butadiene 


argon. 
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Figure 49a. Comparison of the amount of soot formed at a 

reaction time of 1.0 ms for the mixtures 0.311% 
benzene-0. 54% 1,3-butadiene-argon, 0.311% benzene- 
argon, 0.311% benzene-1 .09% acetylene-argon, and 
0.54% 1,3-butadiene-argon. 
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Figure 50. Soot induction times for the mixtures 0.54% 
1,3-butadiene-argon, 0.311% benzene-argon, 
0.311% benzene-0. 54% 1,3-butadiene-argon and 
0.311% benzene- 1 .09% acetylene-argon. 
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Figure 51: Soot formation rates for the mixtures 0.54% 

1,3-butadiene-argon, 0.311% benzene-argon, 
0.311% benzene-0 .54% 1,3-butadiene-argon 
and 0.311% benzene- 1 . 09% acetylene-argon. 
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DISCUSSION OF TEST RESULTS 
1. Effect of Oxygen 

The effects produced by the addition of relatively small amounts of 

molecular oxygen on soot formation from allene (Figs. 9-12), 1,3- 

butadiene (Figs. 13-16), vinylacetylene (Figs. 18-22) and chlorobenzene 

(Figs. 31-35) have been investigated. Analysis of these results 

indicates that the effect of oxygen is principally similar to the effect 

1 2 

observed for acetylene and toluene * . Soot promotion at lower 
temperatures and its suppression at higher ones can be explained as 
previously by interaction between oxidizing agents and active 
intermediates’ at lower temperatures the addition of oxygen enhances 
the rate of fuel pyrolysis while at higher temperatures soot precursors 
are destroyed by oxidation. 

There are, however, a number of distinct features. First, in the 
case of vinylacetylene, the addition of oxygen results in an increase of 
the maximum amount of soot formed, which is particularly noticeable at 
longer observation times (Fig. 20). For all other hydrocarbons studied, 
the addition of oxygen always reduced the peak amount of soot at a given 
reaction time. It is also interesting to note that the induction time 
for soot appearance has a significantly different temperature dependence 
in the case of vinylacetylene oxidation compared to that of its 
pyrolysis (Fig. 21). These dependences for other studied hydrocarbons 
are usually very similar (see e.g Figs. 11, 15 and 34). 

A second distinct feature is the much smaller effect of oxygen on 
soot formation from chlorobenzene compared to that from benzene (Figs, 
28,32,33). It can be explained, as will be discussed in more detail 
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later, by the faster pyrolysis initiation reaction of chlorobenzene 
compared to benzene. For a faster pyrolysis the influence of 
competitive reactions, in this case various oxidation reactions, should 
be smaller. 


2. Effect of Fuel Structure 
a) Pyrolysis of individual hydrocarbon s 
Figures 5-8, 17, and 23-30 present the results and their 

intercomparisons for additional hydrocarbons studied during the reported 
period. Figure 5 depicts the soot yields obtained in pyrolysis of 
ethylene and Figs. 6-8 compare these results with those of acetylene. 
The lower sooting tendency of ethylene compared to acetylene is in 
accord with the results of Fussey and co-workers^ ’ ^ . Fussey et al.^ 
explained the difference "by the variation in the time required to 
produce polyu-'^tvlenes from the reactant", assuming polyacetylenes to be 
the major soot precursors in that "the various hydrocarbons are likely 
to have a similar route to carbon from these species onwards". 
Polyacetylenes were not found to be important intermediates for soot 
formation in our modeling study, which will be discussed later in the 
report. Nevertheless, the route to soot identified by the model would 
be initiated faster in the pyrolysis of acetylene than in the pyrolysis 
of ethylene, which explains our and Fussey et al.'s experimental 
results . 

Following similar arguments, the results obtained with 1-butene 
(Figs. 23-26), which are very close to those with 1,3-butadiene, can be 
rationalized. What is surprising, however, is the relatively low 
sooting tendency of vinylacetylene compared to allene (Figs. 17, 2**) . 




62 


12—17 

Various theoretical considerations indicate vinylacetylene to be as 

effective a soot precursor as allene. Our results, however, clearly 

show that soot formation from vinylacetylene is comparable to that from 

1,3-butadiene and 1-butene rather than to that from allene. 

Figures 27-30 present the results obtained in pyrolysis of a 

chlorobenzene mixture, where they are compared to the results obtained 

at similar conditions in the benzene mixture. Examining Fig. 28, one 

observes that the soot-vield bell is shifted to lower temperatures by 

about 150 K compared to benzene. The shift and its value can be 

explained by the difference between the cleavage of C-H and C-Cl bonds 

in benzene and chlorobenzene molecules, respectively. 

The position of the soot-yield maximum in pyrolysis of aromatic 

5 18 

hydrocarbons is determined ’ by the ratio 

r = k f /k p [A) 0 , (1) 

wheie k f is the rate constant of ring fragmentation, [A]q is the initial 
hydrocarbon concentration, and k is the rate constant of the 

p 

polymerization process forming soot. Assuming the maximum appears at 
the same value of r, one can write 

kf/k p tA] 0 - k f /k p [A ] Q , (2) 

where the left-hand side of Eq. (2) is ratio (1) for benzene and the 
right-hand side is that for chlorobenzene. For the conditions used in 
this study [A]q = [ A ' ] q (see Table I) and assuming k p = k p , Eq. (2) can 
be rewritten as 

k f = kj. (3) 

3 19 

Considering the energetics * of the elementary reactions which may 
constitute the fragmentation path for benzene molecule, the following 


sequence 
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0H •» 0 + H M l800 =464 kJ/mol; ASj g()0 =153 J/mol-K, (Rl ) 

0 -*> HChC-CH=CH-CH=CH AHj g00 =247 kJ/mol; AS° g00 = 65 J/mol-K, (R2) 

where 0 denotes phenyl radical, is suggested for the experimental 
conditions used in this study. The rate limiting step of this sequence, 
particularly during the initial phase of pyrolysis, is reaction (Rl). 
Hence, the rate constant of this reaction, k^, should be considered for 
the left-hand side of Eq. (3). By analogy, the right-hand side of Eq. 
(3) is s'-t equal to the rate constant of the following reaction 

0C1 -» 0 + Cl AhJ 800 =365 kJ/mol ; AS° 80() =123 J/mol-K 20 ’ 21 . (R3) 

Assuming the Arrhenius form for the rate coefficient expressions and 
assuming the ratio of the rate constants for the reverse directions of 

I/O 

reactions (Rl) and (R3) k_/k^ = (m'/m) where |J and (j ' are the 

reduced masses for the corresponding recombinations, Eq. (3) can be 
rewritten as 

R£n(p7p) + (AS-AS') = AH/T - AH'/T' (4) 

mm 

where and are the temperatures of soot yield maxima in mixtures M 

and P, respectively. For a reaction time of 1 ms, is approximately 

equal to 1900 K (Fig. 28). Using the above thermodynamic data for 1800 

K, the approximate midpoint of the temperature interval of interest, we 

obtain from Eq. (4) T' = 1772 K, so that the temperature shift, T -T', 

m mm 

equals approximately 130 K. This is in good agreement with the 
corresponding experimental value. 

The above discussion is in disagreement with the conclusion of 
22 

Singh and Kern that "the phenyl radical is not the primary product of 
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benzene decomposition... (and] ... fragmentation of benzene (presumably 

by an elementary reaction] is the dominant pathway of the decomposition 

in the temperature range of 1400-2200 K" . Their conclusion is based on 

23 

the observation, previously noted by Smith and Johnson , that a larger 
concentration of phenyl radicals is observed in the pyrolysis of 
chlorobenzene than in the pyrolysis of benzene at similar conditions. 
This result can be explained, however, by a simple kinetic argument of 
consecutive reactions. Indeed, in a two-step sequence 

k l k 2 

0H ■* 0 ■+ products , 

the concentration of the intermediate, 0, increases with the ratio 
k^/k^. Based on the previous discussion, for chlorobenzene this ratio 
at 1800 K is approximately 8 times larger than a similar ratio for 
benzene, which should result in a higher concentration of phenyl 
radicals for the former case. 

b) Pyrolysis of hydrocarbon mixtures 
Figures 36-51 present the results obtained in pyrolysis of binary 
hydrocarbon mixtures. Figure 36 presents soot yields obtained in a 
acetylene-hydrogen mixture and Fig. 37 depicts the comparison of these 
results with those obtained at similar conditions in the acetylene 
mixture. As can be seen in this figure, hydrogen strongly suppresses 
soot formation from acetylene. Our chemical kinetic model, which will 
be discussed later, predicts such behavior. The reason for the soot 
suppression, as follows from the model, is the increased rate of reverse 
reactions of hydrogen atom. 

The effecc of addition of acetylene on soot formation from benzene 
and chlorobenzene is examined in Figs. 38-47. It is instructive to 
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examine the effect of a hydrocarbon additive not only by analyzing 
fractional soot yields, as has been the usual practice in our 
works* as well as others^’^’^, but also by considering the 
actual amounts of soot formed. For this purpose, we kept the molecular 
concentration of the main reactant, e.g., benzene, constant in Mixtures 
M, P, S, T, U, V and W (see Table I) rather than the total concentration 
of carbon atoms. While both approaches provide the same information for 
individual fuels, the information is complementary for fuel mixtures 
(cf. Figs. 39 and 39a, Figs. 45 and 45a, Figs. 49 and 49a). Thus, in 
the case of chlorobenzene, analysis of Fig. 39a indicates that the 
addition of acetylene practically does not affect the amount of soot 
formed at low temperatures, while Fig. 39 shows that the fractional 
conversion to soot is not changed at high temperatures. The 1 ^ter 
result is easily understood: at high temperatures the aromatic rings 
are entirely fragmented with acetylene being the main product. The 
only thing that matters is the total amount of carbon atoms (or 
acetylene concentration) present. 

The low-temperature observation in Fig. 39a is not as simple as 
it may first appear - no influence of acetylene addition on soot 
formation from chlorobenzene. Indeed, our computational results, which 
are discussed later in the report, indicate that acetylene addition to 
aromatic radicals is one of the most important reactions in the poly- 
merization growth of aromatic hydrocarbons. Why then does the addition 
of acetylene have no effect? While the answer to that question is 
presently under investigation in our laboratory, a preliminary expla- 
nation can follow from analysis of the rest of the experimental data. 

Inspection of Figs. 40 and 41 indicates the addition of acetylene 
to chlorobenzene increases induction times for soot appearance (Fig. 



66 


40) and also increases the rates of soot production (Fig. 41). In other 

words, acetylene inhibits soot formation during the induction period 

and enhances it after that. The inhibition can be explained by the 

reaction between phenyl radical and acetylene: being fast and 

producing a relatively stable species (e.g. phenylacetylene) , this 

reaction competes with phenyl fragmentation. Since the radical HChC- 
* 

CH=CH-CH=CH, formed in the fragmentation, and immediate products of 
its decomposition., species, are presumably more efficient 

"building blocks" than acetylene^, the removal of phenyl radical from 
the fragmentation route may explain the observed phenomenon. 

That the acetylene reactions are in a fine, competing balance is 
demonstrated by the results of acetylene addition to benzene (Figs. 42- 
47). At low acetylene concentration, there is no significant change in 
the amount of soot formed (Fig. 45a). However, as the concentration of 
acetylene is increased, the enhancement of soot production is well- 
pronounced. The results obtained with a benzene-1 ,3-butadiene mixture 
(Figs. 48-51) further support the hypothesis that C^H x species may be 
more efficient building blocks than acetylene: the addition of 1,3- 

butadiene increases both the fractional soot yields (Fig. 49) and the 
absolute amounts of soot (Fig. 49a) formed from benzene at high tem- 
peratures, indicating a synergistic effect. 

It is interesting to note that the increase in induct on time with 
the addition of acetylene is more pronounced in the case of chloro- 
benzene (Fig. 40) than in the comparable case of benzene (Fig, 46). The 
reason for this is the larger rate of production of phenyl radicals 


in the former case. 


* 
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DETAILED KINETIC MODELING 
1. Introduction 

Soot formation is an important and persistently investigated aspect 
15 16 25-27 

of combustion ’ ’ Nonetheless, the elementary chemical reac- 

tions leading to soot formation are still unknown. There is growing 
evidence that the reactions leading to soot formation involve only 
carbon- and hydrogen-containing species, as oxygen is bound in rich 
flames as CO and H^O. For this reason a large number of studies of soot 
formation in hydrocarbon pyrolysis have been undertaken using flow 
reactors and shock tubes. While these experiments have not provided 
detailed composition profiles, they cover wide ranges of conditions and 
are more easily subjected to interpretive modeling studies at this time. 

The objective of the present modeling was to discover the main 
chemical reaction pathways to soot by experimenting with detailed 

kinetic models of soot formation under the conditions used in shock-tube 

24 

pyrolysis of acetylene . The approach taken was to develop a mechanism 
composed of conventional elementary reactions and to compare the 
predicted time scale of soot formation and absolute values of soot 
yields with their experimental counterparts. Nearly all of the 
equilibrium constants and rate coefficients had to be estimated. A 
method to account for the infinite growth of soot mass also had to be 
developed . 

After the reaction paths implied by our assumed thermochemistry and 

rate constants had been identified by study of the simulations, we 

28 

became aware of a paper by Bockhorn et al. in which astonishingly 
similar conclusions had been drawn from consideration of experimental 
composition profiles in sooting flames. 
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This part of our work was performed in collaboration" 5 with 
Professor William C, Gardiner, Jr., of the University of Texas and Dr. 
Stephen E. Stein of the National Bureau of Standards. 

2. Computer Model 

The reaction aechanism has three logical components: a set of 
reactions describing acetylene pyrolysis, a set of reactions describing 
the formation of larger molecules and radicals and eventually small 
aromatic molecules and radicals, and a description cf the further growth 
of aromatic rings. A total of approximately 600 elementary reversible 
reactions of 180 species were considered during the course of this 
modeling study. 

The chemical species used in the mechanism are listed in Table II 
with both ♦'heir structures and assigned computer names. Atonic hydrogen 
and molecular hydrogen, whose computer names (H and H2) are identical 
with those used conventionally, are not listed in Table II. The names 
were chosen for two principal reasons: first, complete IUPAC names 
were too long and unwieldy; and second, each nan"' was required to 
provide a reasonable description of the species' structure. Our nomen- 
clature system is not an attempt to design a universal method; it was 
developed to address the needs of our specific modeling work ca soot 

formation. The task of developing a completely general nomenclature 

29 

system is nontrivial 

The chemical species are described using 12 characters or less. 
The followirg terms are used to specif} the basic structure of a 
compound : 
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TABLE II. REACTION SPECIES. 
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An: species of type A containing n fused aromatic rings. Species 

with these names include highly condensed polyaromatics such 
as pyrene, A4, and coronone, A7, (see Table II) and the 
intermediates (of Route 1 in Fig. 53) leading to their 
formation, such as phenanthrene , A3, and benzofa] pyrene, A5 . 
Species of type A with four or more rings are both ortho- and 
peri-fused; that is, the aromatic rings share i common faces 
but less than 2i common carbon atoms. 

Bn: species of type B containing n fused aromatic rings. These 

species are involved in the external cyclization route (Route 
2 of Fig. 53). Examples of species with these names are 
benzo[c]phenanthrene , B5 , and dibenzo[c ,d Jphenanthrene , B6 
(see Table II). Species of type B with fewer than seven rings 
are ortho-fused; i.e., the aromatic ritgs share i common sides 
and 2i common carbon atoms. 

Rn: non-aromatic ring containing n carbon atoms. 

CnHm: non-ring structure containing n carbon atoms and m hydrogen 

atoms; neither n nor m is printed if it equals one. 

If a structure contains both aromatic and nonaromatic rings, the 
aromatic rings are specified first. For example, acenaphthalene is 
named as A2R5 (see Table II). 

Ring substituents are listed as suffixes. They are listed in order 
of increasing number of carbon atoms; free-radical substituents come 
last. If two groups have an equal number of carbon atoms, they are 
listed in order of increasing number of hydrogen atoms. Multiple iden- 
tical substituents are denoted by placing a right parenthesis and the 
number of such groups after the substituent name. 
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The names of all free radicals, other than ethynyl radicals, 
contain a special symbol indicating the type of radical. The symbol is 
placed after the structure description. These symbols are: 

13: a vinyl radical with the unpaired electron located on a 

terminal carbon atom. Example: C4H3U (see Table II). The 

only exception to this naming rule is for which was 

named C2H3 in accordance with convention. 

S: a vinyl radical with the unpaired electron located on an 

internal carbon atom. Example: C4H3S (see Table II). 

*: a radical with the unpaired electron located on an aromatic 

ring and adjacent to a C2H substituent. Example: A1C2H* (see 

Table II). 

-: a radical which has the unpaired electron located on an 

aromatic ring and which is not described by Example: 

Al- (see Table II). 

%: a radical with the unpaired electron located on a non-aromatic 

ring. Example: R6% (see Table II). 

If more than one structure would be assigned the same name using 
these tul^c, thr structures are distinguished by letters placed at the 
ends of the names; e.g. the second species name would end with X, the 
third would end with Y and so on. When our initial computer experi- 
mentation established that distinctions between particular isomers were 
not important, the isomers were lumped under one name. 

The reaction set is presented in Table III. The core of the 

acetylene pyrolysis model is based on the mechanism of Tanzawa and 
■an 

Gardiner whose development entailed comparison with all of the 
available data. It was thus calibrated overall against the measured 
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H 

+ 

C8H5S 

■ 

C8H6 



29 

-426.8 

( 86) 

H2 


C8H5S 

m 

H 

4 

C8H6 

-2 

21.8 

< 87) 

C2K2 

+ 

C8H5S 

m 

C2H 

4 

C8H6 

-4 

79. 1 

( 88) 

H 

+ 

Al- 

* 

A1 



21 

-470.8 

( 89) 

A 1 

■ 

ii 

4 

C6H5U 



22 

673.0 

( 90) 

C2H 

+ 

Al- 

« 

A1C2H 



7 

-512.4 

( 91) 

A1C2H 

m 

ii 

4 

C8HSSX 



22 

619.8 

< 92) 

H 

+ 

A1C2H- 

» 

A1C2H 



21 

-470 . 8 

< 93) 

H2 

+ 

A1C2H- 

m 

H 

4 

AIC2H 

-23 

-22. 1 

( 94) 

C2H2 

■f 

A1C2H- 

*a 

C2H 

4 

AIC2H 

-24 

35.2 

( 95) 

C2H2 

+ 

A1C2H- 

m 

H 

4 

A 1 C2H ) 2X 

-25 

-6.5 

( 96) 

H2 

+ 

A 1 - 

B 

H 

4 

A 1 

-23 

-22. 1 

< 97) 

C2K2 

+ 

A 1 - 

m 

C2H 

4 

A) 

-24 

35.2 

< 98) 

C2H2 

4 

A 1 - 

m 

H 

4 

AIC2H 

-25 

-6.5 

< 99) 

C4H2 

+ 

A 1 - 

m 

C4H 

4 

A 1 

-24 

37.2 

<100) 

C4H2 

+ 

A 1 - 

-*» 

C2H 

4 

AIC2H 

26 

57.0 

<101 ) 

C4H2 

4 

Al- 

o 

A1C4H2S 



27 

-181.6 

<102) 

H 

+ 

A1C4H2S 

» 

A1C4H3 



29 

-426.8 

<103) 

A1C4H3 

“ 

C4H3U 

4 

Al- 



-7 

455.7 

<104) 

H2 

+ 

A1C4H2S 

■ 

H 

4 

A1C4H3 

-2 

21.8 

<105) 

C2H2 

+ 

A1C4H2S 

m 

C2H 

4 

AIC4H3 

-4 

79. 1 

<106) 

C2H2 

+ 

A1C4H2S 

m 

H 

4 

AIC6H3 

1 1 

17.7 

<107) 

C4H2 

4 

A1C4H2S 

* 

A1C8H4S 



.3 

-219.0 

<108) 

A1C8H4S 

at 

K 

4 

A2C2H >2X 



20 

44.2 

<109) 

H 

4 

A1C8H4S 

■ 

A 1 C8H5 



29 

-426.8 

< 1 10) 

H2 

4 

A1C8H4S 

m 

K 

4 

A1C8H5 

-2 

21.6 

<111) 

C2H2 

4 

A 1 C8H4S 


C2H 

4 

A1C8H5 

-4 

79. 1 

< 1 12) 

C2H2 

4 

A1C8H4S 

tn 

H 

4 

AICI0H5 

1 1 

17.7 

<113) 

C2H2 

4 

Al- 

■ 

A 1 C2H2U 



27 

-170.0 

<1 14) 

A1C2H2U 

B 

ii 

4 

A1C2H 



-5 

163.5 

<115) 

K 

4 

A1C2H2U 

M 

A1C2H3 



29 

-460.3 

< 116) 

H2 

4 

A1C2H2U 

a 

H 

4 

A1C2H3 

-2 

-11.6 

< 117) 

C2H 

4 

A1C2H2U 


A1C4H3 



7 

-50 1 . 9 

<118) 

C2H2 

4 

A1C2H2U 

m 

A1C4H4U 



13 

-22 1 . 1 

<119) 

C4H2 

4 

A 1 C2H2U 

m 

A1C6H4S 



13 

-232.7 

020) 

A1C4H4U 

s 

H 

4 

A2 



20 

-8,9 

<121 ) 

A1C4H4U 

K 

H 

4 

A1C4H3 



-5 

225. 1 

022) 

H 

4 

A1C4H4U 

ac 

A1C4H5 



29 

-460.3 

023) 

K2 

4 

A1C4H4U 

A 

H 

4 

A1C4H5 

-2 

-11.6 

024) 

C2H2 

4 

AIC4H4U 

m 

C2H 

4 

A1C4H5 

-4 

45.6 

< 125) 

C2H2 

4 

AIC4H4U 

e 

H 

4 

AIC6H5 

1 1 

4.0 

< 126) 

A1C6H4S 

s 

H 

4 

A2C2H 



20 

24.5 

< 127) 

A1C6H45 

B 

H 

4 

A1C6H3 



-5 

238.8 

< 128) 

K 

4 

A1C6H4S 

« 

A 1 C6H5 



29 

-426.8 

< 129) 

K2 

4 

A1C6H4S 

■ 

H 

4 

A1C6K5 

-2 

21 .8 
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(130) C2H2 

+ A1C6H4S 

- C2H 

+ A1C6H5 

-4 

79. 1 

(131) C2H2 

+ A1C6H4S 

- H 

+ A1C8H5 

I 1 

17.7 

(132) C2K2 

+ A1C4H2S 

■ AIC6H4U 


13 

-135.1 

(133) A1C6H4U 

■ H 

+ A2C2H 


20 

-61.5 

(134) A1C6H4U 

- H 

+ A1C6H3 


-5 

152.8 

(135) H 

+ A1C6H4U 

• A1C6H5 


29 

-512.8 

(136) H2 

+ A1C6H4U 

■ H 

+ A1C6H5 

-2 

-64.2 

(137) C2H2 

+ A1C6H4U 

- C2H 

+ A1C6H5 

-4 

-6.9 

(136) C2H2 

+ A1C6H4U 

- H 

+ A1C8H5 

1 1 

-66.3 

(139) H 

+ A2- 

« A2 


21 

-470.8 

(140) H2 

+ A 2- 

■ H 

+ A2 

-23 

-22. 1 

(141) C2H 

+ A2- 

* A2C2H 


7 

-512.4 

(142) C2H 2 

+ A2- 

• C2H 

+ A2 

-24 

35.2 

(143) A1C2H 

- H 

+ A 1C2H« 


-21 

470.8 

(144) H 

+ A1C2H 

- H2 

+ A1C2H« 

23 

22. 1 

(145) C2H 

+ AIC2H 

- C2H2 

+ A1C2H* 

24 

-35.2 

(146) n)C2H- 

■ C8H5U 



-30 

182.5 

(147) H 

+ C8H5U 

« C8H6 


29 

-460.3 

(146) H2 

+ C8H5U 

* H 

+ C8H6 

-2 

-11.6 

(149) C2H2 

+ C8H5U 

*= C2H 

+ C8H6 

-4 

45.6 

(150) C2H2 

+ C8HSG 

“ H 

+ CI0H6 

1 1 

-15.8 

(151) C8H5U 

« H 

+ C8H4 


-E 

205.3 

(152) C8H5U 

- C2H2 

+ C6H3U 


-13 

221.1 

(153) H 

+ C6H3U 

* C6H4 


29 

-460.3 

(154) H 2 

+ C6M3U 

«* H 

+ C6H4 

-2 

-11.6 

(155) C2H2 

+ C6H3U 

■ C2H 

+ C6H4 

-4 

45.6 

(156) C2H2 

+ C6H3U 

- H 

+ C8H4 

1 1 

-15.8 

(157) C6H3U 

■= C2H 

+ C4H2 


-8 

214.2 

(158) C2H2 

+ A1C2H* 

“ A1C2HC2H2U 


27 

-170.0 

(159) C2H2 

+ A 1 C2H» 

= H 

+ A 1 C2H )2 

-25 

-6.5 

(160) A 1C2HC2H2U 

« H 

+ A1C2H)2 


-5 

163.5 

(161) A IC2HC2H2L' 

*= A2-X 



30 

-223.5 

(162) H 

+ A 1C2HC2H2U 

= A1C2HC2H3 


29 

-460.3 

(163) H2 

+ A 1 C2HC2H2U 

* H 

+ A1C2HC2H3 

-2 

-11.6 

(164) C2H2 

+ A1C2HC2H2U 

= C2H 

+ A IC2HC2K3 

-4 

45.6 

( i65) H 

+ A2-X 

*= A 2 


21 

-470 . 8 

(166) H2 

+ A2-X 

- H 

+ A2 

-23 

-22. 1 

(167) C2H2 

+ A2-X 

= C2K 

+ A2 

-24 

35.2 

(168) C2H 

+ A2-X 

- A2C2H 


7 

-512.4 

(169) C4H2 

+ A1C2H« 

- AIC2HC4H2S 


27 

-181.6 

(170) A1C2HC4H2S 

*= H 

+ A1C2HC4H 


-5 

177.2 

(171) A1C2HC4H2S 

*= A2C2H- 



30 

-190.1 

(172) H 

+ A1C2HC4H2S 

-= A1C2HC4H3 


29 

-426 . 8 

(173) H2 

+ A1C2HC4H2S 

■ H 

+ A1C2HC4H3 

-2 

21 .8 

(174) C2H2 

+ A 1C2HC4H2S 

- C2H 

+ AIC2HC4H3 

-4 

79. 1 

(175) H 

+ A2C2H-X 

= A2C2H 


21 

-470.8 

(176) H2 

+ A2C2H*X 

«= H 

+ A2C2H 

-23 

-22. 1 

077) C2H2 

+ A2C2H»X 

*= C2H 

+ A2C2H 

-24 

35.2 

(178) C2H2 

+ A2C2K*X 

= H 

+ A2R5C2H 

27 

-134.2 

(179) H 

+ A 1 

« R6X 


19 

-90.7 

(160) R6X 

= C6H7U 



-20 

78.4 

(181) C2H2 

+ A2-X 

- H 

+ A2C2H 

-25 

-6.5 

(182) C2H2 

+ A2- 

■ H 

+ A2C2H 

-25 

-6.5 

(183) C2H2 

+ A2C2H- 

■ C2H 

+ A2C2H 

-24 

35.2 

084) H 

+ A2C2H* 

*= A2C2H 


21 

-470.8 

0 85) K2 

+ A2C2H" 

= H 

♦ A2C2H 

-23 

-22.1 

086) C2K2 

+ A2C2H* 

= H 

+ A2C2HJ2 

-25 

-6.5 

087) C2H2 

+ A2C2H- 

■ A2C2HC2H2U 


27 

-170.0 

(188) A2C2HC2H2U 

= H 

+ A2C2H ) 2 


-S 

163.5 

089) A2C2HC2H2U 

= A3- 



30 

-232.7 

0 90) H 

+ A3- 

- A3 


21 

-470. S 

091) H2 

+ A3- 

*= H 

+ A3 

-23 

-22. 1 

092) C2H2 

+ A3- 

- C2H 

+ A3 

-24 

35.2 

0 93) C2H2 

+ A3- 

- H 

+ A3C2H 

-25 

-6.5 

094) C2H2 

+ A3- 

*= A3C2H2U 


27 

-170.0 
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(195) 

A3C2H2U 

- H 

+ A4 


20 

-50.6 

( 196) 

A2- 

- A1C4H3- 



-30 

229. 1 

( 197) 

H 

+ A1C4H3- 

- A1C4H3 


21 

-465 . 8 

( 196) 

H2 

+ A1C4H3- 

- H 

+ A1C4H3 

-23 

-17.2 

< 199) 

C2H2 

+ A1C4H3- 

- C2H 

+ A1C4H3 

-24 

40. 1 

(200) 

C2H2 

+ A1C4H3- 

« H 

+ AIC2HC4K3 

-25 

-1.6 

(201 > 

A1C4H3- 

- C10H7U 



-30 

125.9 

(202) 

C10H7U 

- K 

+ CI0H6 


-5 

205.3 

(203) 

C10H7U 

■ C2H2 

+ C8H5U 


-13 

221.1 

(204) 

H 

+ A4- 

» A4 


21 

-470.6 

(205) 

H2 

+ A4- 

- H 

+ A4 

-23 

-22. 1 

(206) 

C2H2 

+ A4- 

- C2H 

+ A4 

-24 

35.2 

(207) 

C2H2 

+ A«- 

- H 

+ A4C2K 

-25 

-6.5 

(208) 

H 

♦ A*'C2H» 

- A4C2H 


21 

-470.8 

(209) 

H2 

+ A4C2H" 

- H 

+ A4C2H 

-23 

-22. 1 

(2)0) 

C2H2 

♦ A4C2K- 

» C2H 

+ A4C2H 

-24 

35.2 

(21 1 ) 

C2H2 

+ A4C2H* 

- H 

+ A4C2H)2 

-25 

-6.5 

(212) 

C2H2 

+ A4C2H" 

« A4C2HC2H2U 


27 

-170.0 

(213) 

A4C2HC2H.2U 

- H 

+ A4C2H)2 


-5 

163.5 

(214) 

A4C2HC2H2U 

- A5- 



30 

-241 .9 

(215) 

H 

+ A5- 

“ A5 


21 

-470.8 

(216) 

H2 

+ A5- 

- H 

+ A5 

-23 

-22. 1 

(217) 

C2H2 

+ A5- 

- C2H 

+ A5 

-24 

35.2 

(218) 

C2H2 

+ A5- 

- H 

+ A5C2H 

-25 

-6.5 

(21* ) 

C2H2 

+ A5- 

** H 

4 A6 

27 

-220.6 

(220) 

H 

+ A6- 

■ A6 


21 

-470.8 

(221 ) 

H2 

+ A6- 

- H 

+ A6 

-23 

-22. 1 

(222) 

C2H2 

+ A6- 

- H 

+ A6C2H 

-25 

-6.5 

(223) 

C2H2 

+ A6- 

- C2H 

+ A6 

-24 

35.2 

(224) 

C2H2 

+ A6- 

* H 

+ A7 

27 

-220 . 6 

(225) 

H2 

+ A3C2H- 

» H 

+ A3C2H 

-23 

-22. 1 

(226) 

H 

+ A3C2H* 

« A3C2H 


21 

-470.8 

(227) 

C2H2 

+ A3C2H- 

■ C2H 

+ A3C2K 

-24 

35.2 

(228) 

C2H2 

+ A3C2H* 

«■ H 

+ A3C2H>2 

-25 

-6.5 

(229) 

C2H2 

+ A3C2H* 

*>. A3C2HC2H2U 


27 

-170.0 

(230) 

A3C2HC2H2U 

- H 

+ A3C2H ) 2 


-5 

163.5 

(231 ) 

A3C2HC2H2U 

■ B4- 



30 

-203.4 

(232) 

H2 

+ B4- 

- H 

4 B4 

-23 

-22. 1 

(233) 

H 

+ B4- 

“ B4 


21 

-470.8 

(234) 

C2H2 

+ B4- 

- C2H 

4 B4 

-24 

35.2 

(235) 

C2H2 

+ B4- 

- H 

4 B4C2H 

-25 

-6.5 

(236) 

K2 

+ B4C2H* 

« H 

4 B4C2H 

-23 

-22. 1 

(237) 

H 

+ B4C2H* 

■ B4C2H 


21 

-470.6 

(238) 

C2H2 

+ B4C2H* 

- C2H 

4 B4C2H 

-24 

35.2 

(239) 

C2H2 

+ B4C2H* 

- H 

+ B4C2H >2 

-25 

47.9 

(240) 

C2H2 

+ B4C2H« 

« B4C2HC2H2U 


27 

-115.6 

(241 ) 

B4C2HC2H2U 

*= H 

+ B4C2K ) 2 


-5 

163.5 

(242) 

B4C2HC2H2U 

= B5- 



30 

-274 . 5 

(243) 

B5- 

- H 

+ A6 


27 

-83.0 

(244) 

H 

+ B5- 

■ B5 


21 

-428.9 

(245) 

H2 

+ B5- 

* H 

+ B5 

-23 

19.7 

(246) 

C2H2 

+ B5- 

•* C2H 

+ B5 

-24 

77.0 

(247) 

A2C2H* 

» A 1 C2H» 



31 


(248) 

A4C2H- 

- A3- 



31 


(249) 

K 

4 A7- 

- A7 


21 

-470.8 

(250) 

K2 

+ A7- 

* H 

4 A7 

-23 

-22. 1 

(251) 

C2H2 

+ A7- 

b C2H 

4 A7 

-24 

35.2 

(252) 

C2H2 

+ A7- 

- H 

4 A7C2H 

-25 

-6.5 

(253) 

H 

+ A7C2K« 

» A7C2H 


2) 

-470.8 

(254) 

H2 

+ A7C2H" 

« H 

4 A7C2H 

-23 

-22. 1 

(255) 

C2H2 

+ A7C2H- 

* C2H 

4 A7C2H 

-24 

35.2 

(256) 

C2H2 

+ A7C2H- 

B H 

4 A7C2HJ2 

-25 

-6.5 

(257) 

C2H2 

+ A7C2H* 

B A7C2HC2H2U 


27 

-170.0 

(258) 

A7C2HC2H2U 

« H 

+ A7C2H )2 


— D 

163.5 

(259) 

A7C2HC2H2U 

* A8- 



30 

-241.9 
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(260) H 

+ A8- 

- A8 

(261) H 2 

+ A8- 

- H 

<262) C2H2 

+ A8- 

- C2H 

<263) C2H2 

+ A8- 

- H 

<264) C2H2 

+ A8- 

- H 

<265) H 

+ A9- 

- A9 

<266) H2 

+ A9- 

- H 

<267) C2H2 

+ A9- 

- H 

<266) C2H2 

+ A9- 

- C2H 

<269) C2H2 

+ A9- 

- H 

(270) C2H2 

+ A2-X 

- H 

<271) H 

+ A2R5- 

■ A2R5 

<272) C2H2 

+ A2R5- 

- C2H 

<273) H2 

A2R5- 

- H 

<274) C2H2 

+ A2R5- 

- H 

<275) H 

+ A2R5C2H 

- H2 

<276) H 

+ A2R5C2H" 

- A2R5C2H 

<277) C2H2 

+ A2R5C2H* 

- C?H 

<276) C2H2 

+ A2R5C2H" 

* H 

<279) C2H2 

+ A2R5C2H" 

- A2R5C2HC2H2U 

<260) A2R5C2HC2H2U 

- H 

+ A2R5C2H)2 

<->31 ) A2R5C2HC2H2U 

- A3RS- 


<282) H 

+ A3R5- 

- A3R5 

<283) H2 

+ A3RS- 

- H 

<284) C2H2 

+ A3R5- 

- C2H 

<265) C2H2 

+ A3R5- 

- A3R5C2H2U 

<286) A3R5C2H2U 

- H 

+ A3R5C2H 

<287) A3R5C2H2U 

■ H 

+ A4 

<288) C2H2 

+ C6H5U 

■ H 

<289) C2H 

+ C8H5SX 

* C10H6 

<290) C2H2 

+ C8H5SX 

■ H 

<291 ) C4H 

+ C6R 

- C10K2 

<292) C6H 

+ C6H 

- CI2H2 

<293) C6H 

+ C6H2 

■ C12H3S 

<294) C4H 

+ C6H3S 

- C10H4 

<295) C2H3 

+ C6H3S 

■ C8H6' 

<296) C4H3S 

+ C6H3S 

* C10H6 

<297) C6H3S 

+ C6H3S 

- C12H6 

<298) C6H2 

+ C6H7S 

« C2H 

<299) C6H2 

+ C6H3S 

« C12H5S 

<300) C2H 

+ C8H3S 

» C10H4 

<30 1 ) C2K3 

+ C8H3S 

- C10H6 

<302) C4H3S 

+ C8H3S 

- C12H6 

003) C6H3S 

+ C8H3S 

- C14H6 

004) C2H2 

+ C8H3S 

- H 

005) C4H2 

+ C8H3S 

- C2H 

006) C4H2 

+ C8H3S 

« C12H5S 

007) C6H2 

+ C8H3S 

“ C4H 

008) C6H2 

+ C8H3S 

« C2H 

009) C6H2 

+ C8H35 

- C14H5S 

<310) H 

+ C10H3S 

' CI0H4 

Oil) H2 

+ CI0H3S 

- H 

012) C2H 

+ C10H3S 

* C12H4 

013) C2H2 

+ C10H3S 

- C2H 

014) C2H2 

+ CI0H3S 

- H 

015) H 

+ C12H2 

« C12H3S 

016) C2H 

+ C12H3S 

- C14H4 

017) K 

+ C12H3S 

- CI2H4 
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-4 
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C6H2 

+ A 1 - 
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+ A 1 

-24 

37.2 

(534) 

C6H2 

+ A 1 - 

- C4H 

+ A1C2H 

26 

57.0 

(535' 
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+ AIC4H2S 

- C4H 

+ A1C4H3 

-4 

81.2 
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C4H2 

+ A1C4H2S 

- C2H 

+ A1C6H3 

12 

81.2 

('•37) 

A1C6H3 

- C2H 

+ AIC4H2S 


-7 

468.2 

(538) 

C6H2 

+ AIC4H2S 

- C6H 

+ A1C4H3 

-4 

81.2 
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C6H2 
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12 

81 .2 
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AIC8H5 
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-7 
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= C4H 

+ A1C8H5 

-4 

81.2 

(542) 

C4H2 

+ A1C8H4S 

- C2H 

+ A1C10K5 

12 

81.2 

(543) 

C4H2 

+ AIC8H4S 

- A1CI2H6S 


13 

-2 >9.0 

(544) 
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- C6H 
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-4 

81.2 
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- C4H 
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12 

81 .2 
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C6H4 

+ A 1C8H4S 

= C6H3S 

4 AIC8H5 

-15 

-0.0 

(547) 

C4H4 

+ Al- 

= AIC4H4S 


27 

-253.9 

(548) 

A1C4H4S 
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+ H 

4 A2- 

18 

495.3 

(549) 

H 

+ A 1C4H4S 
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29 

-426.8 

(550) 

A1C4H5 
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+ A1C2H2U 


-7 
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= C2H 

4 AIC4H5 

-4 

79. 1 
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= H 

4 AIC6H5 

1 1 

37.5 

'554) 
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-4 
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= C2H 
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12 

100.9 
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33 

-43.9 

(557) 

C*H4X 

+ A 1 - 
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-0.0 

(571 ) 
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+ A IC6H4SX 
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-0.0 
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4 A 1 C8H5 

12 
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(573) 
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12 

33.5 

(574) 
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= C4H 

4 AIC2H3 

-4 

47.7 

(575) 

C2H 
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(615) 
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H 
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23 
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24 
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A2C2K 

n H 
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-21 
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H 

+ A2C2HJ2X 

- C2H2 

4 

A2C2H- 

25 

6.5 
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chain and non-chain reaction rates. Their mechanism was subsequently 

extended to higher temperatures and also found to be suitable for 

31 

describing ethylene pyrolysis . In our work we distinguish between the 

isomers of the radical, namely H-C=C-C=CH 2 and H-C=C-CH=CH; thus, 

the product of the reaction + ^2^2 = ^4^3 + ^ WaS asslunet ^ t0 t * le 

32 

latter, in line with the theoretical study of Dewar and co-workers 
The rest of the mechanism was generated according to physical organic 
chemistry principles to comprise the likely radical and atom reactions 
(no migrations o i hydrogen atoms or complex isomerizations were 
included) that could lead eventually to cyclization and growth of 
aromatics . 

The rate coefficients for reactions 1-18 were used exactly as given 

in Table II of Ref. 24. Rate coefficients for reactions 19, 20 and 21 

were chosen to equal those of reactions 5, 9 and 8, respectively. For the 

purpose of assigning rate coefficient expressions the reactions were 

grouped into classes as shown in Table IV. The rate coefficient 

expression chosen for a prototype reaction was assigned to the rest of 

the class without change. We assume thereby that a decrease in reaction 

rate due to increasing mass is offset by growing reaction-path 

degeneracy. To assign rate coefficients for classes where no 

experimental data are available, the viewpoint was adopted that 

reasonable upper limits based on collision theory would be preferable to 

transition state theory estimates based on assumed transition structure 

13 

properties. The values of the rate coefficients were limited to 2x10 

3-1 14 

_m mol s for bimolecular reactions of organic species and to 1x10 

3-1 

cm mol s if hydrogen atom is a reactant. The principle of detailed 
balancing was always obeyed. 
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As our modeling indicated rhat it is not important to describe 
some reaction sequences in complete detail, these sequences were 

simplified to keep the reaction mechanism of manageable size. For 

example, reaction class 31, the irreversible fragmentation of aromatic 
radicals, is not a true reaction class. It was used to include ring 
fragmentation effects for large aromatic radicals after detailed models 
of ring fragmentation (including e.g. A2- = A1C4H3-) indicated that 
although ring rupture was an important reaction route for aromatic 
radicals with one ring, it was not important for aromatic radicals with 
more than one ring. Even with this overestimate of fragmentation rates 
for the polycyclic radicals, fragmentation was not a significant 
reaction route for these species. 

Reliable thermodynamic properties were available from the 

literature for only a few of the species of interest. Group additivity 

21 

estimation methods were therefore used to obtain a large fraction of 

necessary thermodynamic data. A detailed discussion of all group values 

used in this work along with illustrations of their applications to 

19 

carbon polymerization reactions is reported elsewhere 

The main test mixture used in the computer experiments was 1.09% 
C 2 H 2 -Ar at ^5 = x 10 5 mol/cm 3 with values from 1500 to 2100 K, 
to simulate Series B of Ref. 24. The pressure range of this series, 
5-7 bar, is consistent with the choice of high-pressure limits for 
unimolecular rate coefficients. The computations were carried out to a 
reaction time of 2 ms using a constant-pressure model. The computed 
temperature changes wer at most a few degrees . The computations were 
done with an IBM 370-3081 computer using the LSODES integrator 3 ^. 


127 


3. Results and Discussion 

The primary objective of our investigation was to discover the 
principal route(s) for formation of soot. Other important aspects, e.g. 
product distributions for various experimental conditions, will be the 
subject of our future studies. Their omission does not affect the 
main conclusions of the present work. 


a) Formation of the first arom . ^ic ring 
Our initial working hypothesis was that a multiplicity of 
cyclization pathways would account for the overall cyclization rate. 
Following this hypothesis, a reaction network was composed that included 
26 different cyclization reactions of species up to C 14 K x' These 
reactions are of two classes, non-radical and radical cyclizations 
(Table IV). An example of the former is (see Ref. 24) 


H x . 

> = c \ 

H - C\_ .C - H 




CH HC^ 


H N • 

> = C \ 

H - C;C - H + H + H 

- C' 


The second cyclization class can be further subdivided into 
a) cyclization of more stable radicals, e.g. 


H - C 5 C 


H H 

= (T 

\ 

C C - H 

H c *. 


H 


H \ 

C = C 

/ \ 

H - C = C - C C-H + H 

V / 

/C " C x 
H % 


TABLE IV. REACTION CLASSES AND RATE COEFFICIENT EXPRESSIONS. 
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b) cyclization of less stable radicals, e.g. 


H YL 

= C 

H \ 

H - C* C - H 

^ // 

/ C ' 

n H 


H 


H - C 


' S C = C / 

/ \ 


\ 


C - H + H 


; c - c x 


H 


and c) radical cyclization via interaction of an unpaired electron with 
a triple bond, e.g. 
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The computational results indicated that for the conditions studied 
cyclization by the latter reaction has 3, 5 and 8 orders of magnitude 
more forward flux than the radical, (a) and (b) , and non-radical 
cyclizations , respectively. The dominant reaction pathway for the 
formation of the first aromatic ring is given in Table V. 

A second important reaction of subclass (c) 
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appeared to proceed rapidly but in the reverse direction, being the 
major fragmentation reaction of the single-ring aromatic hydrocarbons. 

b) Formation of two-ring aromatics 
About 10 essentially different pathways for the formation of two- 
ring aromatic species were considered. As in the first-ring 
cyclization, a single route (Fig. 52) dominated the formation of the 
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TABLE V. FORMATION OF THE FIRST AROMATIC RING. 
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MINOR ROUTE 



i 

H 


Figure 52. Principal reaction pathways for formation of 
two-ring aromatics. For a given reaction step, 
the difference in the lengths of the arrows 
represents the difference in forward and reverse 
reaction rates as computed for the test case 
at T - 1700 K and a reaction time of 0.5 ms. 
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second ring, for the conditions tested, by at least 3 orders of 

magnitude. Also shown in Fig. 52 is the next fastest route * a two-step 

12 

acetylene addition as suggested by Bittner and Howard , which appeared 

to be most important among similar two-step additions such as acetylene 

addition followed by diacetylene addition, diacetylene addition followed 

by acetylene addition, two-step diacetylene addition, or similar 

13 

reactions involving vinyl acetylene . The computational results 
indicated that the two-step addition mechanism is slow because the first 
intermediate adduct, such as phenyl-CH=CH in Fig. 52, is so unstable 
that stabilization via loss of a hydrogen atom is much faster than the 
next addition reaction. 

An important question is why the dominant route is so much faster 
than the rest. The reason is that formation of the more stable 
molecular intermediate phenylacetylene occurs in the dominant route 
(Fig. 52), whereas in the minor routes the molecular species produced 
are relatively unreactive side products and all the reactive 
intermediates are radicals. At first glance one would expect that 
formation of a stable molecule followed by its reactivation via hydrogen 
atom abstraction should be slower than radical-reactant interactions, 
which is true for irreversible kinetics. In the case of consecutive 
reversible reactions the formation of stable intermediates appears to be 
the more decisive factor. 

28 

Recently, Bockhorn et al. suggested a cyclization mechanism very 
similar to ours. Their reaction mechanism was developed to explain 
product distributions obtained in sooting laminar flat flames with a 
variety of hydrocarbons. The agreement in the cyclization steps 
suggested independently and from such different viewpoints as product 


•> 
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distribution in flames in their work and a computer experiment in ours 
is remarkable. 

c) Further growth of fused polycyclic aromatics 

Based on the computational results for two-ring formation, only the 
dominant mechanism of cyclization discussed above was considered for the 
further growth of fused polycyclic aromatics. Two main pathways 
considered in this study are presented in Fig. 53. Although in the 
figure they stop at pyrene (and cyclopenta [cd]pyrene) , both routes were 
continued in the computer experiments to benzo [ghijperylene which was 
followed by formation of coronene. Route 1 in Fig, 53 generally 
dominated the carbon mass flux, although both routes proceed with 
approximately equal rates at the lowest temperatures tested (1500 K). 
At these temperatures, however, the computed soot yields were extremely 
low, as discussed in the next section. Thus, for the growth of fused 
polycyclic aromatics, from coronene to ovalene and further, only Route 1 
of Fig. 53 was considered. 

Routes 1 and 2 in Fig. 53 are similar in that both proceed through 
radical as well as molecular intermediates. The molecular intermediates 
of Route 2 are substituted fused polycyclic aromatics, the unsubstituted 
fused molecules like acenaphthylene, pyrene and coronene having the role 
of relatively unreactive side products. On the other hand, these 
unsubstituted fused aromatics (the cyclopenta group is relatively 
unstable and must be partially lost, in one way or another, along the 
reaction pathway; to indicate this on the diagram of Fig. 53, the group 
is drawn with dashed lines) are included in Route 1 as reactive 
intermediates. What distinguishes these intermediates is their high 
thermodynamic stability, increasing from acenaphthylene to pyrene, to 
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Figure 53. Principal reaction pathways for formation of 
fused polycyclic aromatics. For a given 
reaction step , the difference in tne lengths 
of the arrows represents the difference in 
forward and reverse reaction rates computed 
for the test case at T = 1700 K and a reaction 
time of 0.5 me. 
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coronene and to ovaiene. The computational results indicate that their 
formation reactions are practically irreversible; this is particularly 
true for the formation of acenaphthylene, coronene and ovalene. These 
irreversible reactions have the effect of "pulling" the chain of 
reversible steps, which is crucial for soot formation. 

d) Species concentration profiles 

The time development of the concentrations of selected species at a 
temperature of 1700 K is presented in Figs. 54-58. Figure 54 presents 
the concentrations of H, and H g ^. is the hypothetical hydrogen 
atom concentration which would be in equilibrium with the calculated 
concentrations of H 2 ; that is, (H^] = where K lg is the 
equilibrium constant of reaction 18 (see Table III). The temperature 
dependences of the concentrations of the same species at a reaction time 
of 500 microseconds are shown in Figs. 59-63. 

e) Formation of soot 

For modeling purposes we defined soot, as in our experimental work, 
as a collection of species absorbing light of a specific 
wavelength, or, similarly^, as a collection of species above a certain 
size. In the present modeling the list of these species was chosen to 
begin at coronene. The choice of a larger species would not greatly 
affect the soot yields at longer reaction times. The computed soot 
yield was defined, again as in our experimental work, as the number of 
carbon atoms accumulated in soot divided by the initial number of carbon 
atoms of acetylene. 

To compute soot yields we assumed an infinitely long reaction 
sequence of Route 1 (Fig. 53), all species of this route beginning with 
coienene constituting soot. This is, of course, a simplification: 
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L0G IO (CONCENTRATION), (mol/cm a ) 

12.00 - 11.00 - 10.00 - 9.00 - 0.00 - 7.00 - 0.00 


138 

+ SPECIES H (EQUIL) 
A SPECIES H2 
Q SPECIES H 



1 1 1 j 

10 a . 200 . 300 . UOO . 500 . 

T I ME (m i crosec) 

Figure 54. Time dependence of computed concentrations of 

H, H„ and H for the test case at T = 1700 K. 
I eq 

H is tb.e hypothetical concentration of 
hydrogen atoms which would be in equilibrium 
with the computed concentrations of H,- The 
symbols designate computational results, the 
lines cubic spline fits. 
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X SPECIES C2H3 
+ SPECIES C2H 
A SPECIES H2 



Figure 55. Time dependence of computed concentrations of H, 

C^H and for the test case at T = 1700 K. The 
symbols designate computational results, the lines 
cubic spline fits. 
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Figure 56. Time dependence of computed concentrations of C^H^, 

and for the test case at T = 1700 K. The 

symbols designate computational results, the lines 
cubic spline fits. 
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o SPECIES 
X SPECIES C6H7U 
+ SPECIES C6H5U 
A SPECIES C4H5U 


g o SPECIES C4H3U 

* 



TIME (mi crosec) 

Figure 57. Time dependence of computed concentrations of 

intermediate-sized radicals for the test case at 
T = 1700 K. The structure for each species is 
presented in Table II. The symbols designate 
computational results, the lines cubic spline 
fits . 
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x SPECIES A7 
+ SPECIES m 


A SPECIES A2R5 
O SPECIES R1 



TIME (microsec) 

Figure 58. Time dependence of computed concentrations of 
intermediate- and large-sized aromatic mole- 
cules for the test case at T = 1700 K. The 
structure for each species is presented in 
Table II. The symbols designate computational 
results, the lines cubic spline fits. 
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TEMPERRTURE (K) 


Figure 59. Temperature dependence of computed concentrations 

of H, and h for the test case at a reaction 
Z eq 

time of 0.5 ms. is the hypothetical con- 

centration of hydrogen atoms which would be in 
equilibrium with the computed concentration of 
i^. The symbols designate computational results, 
the lines cubic spline fits. 
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x SPECIES C2H3 
+ SPECIES .C2H 
A SPECIES H2 
O SPECIES H 



TEMPERATURE (K) 


Figure 60. Temperature dependence of computed concen- 
trations of H, H 2 , C^H and for the test 

case at a reaction time of 0.5 ms. The 
symbols designate computational results, the 
lines cubic spline fits. 
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+ SPECIES C6H2 
A SPECIES C4H2 
O SPECIES C2H2 



Figure 61. Temperature dependence of computed concen- 
trations of ^ 2 ^ 2 ’ ^4^? an< * ^6^2 ^° r t ^ ie test 
case at a reaction time of 0.5 ms. The symbols 
designate computational results, the lines 
cubic spline fits. 
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o SPECIES fll- 
x SPECIES C6H7U 
4- SPECIES C6H5U 
A SPECIES C4H5U 
g o SPECIES C4H3U 



Figure 62. Temperature dependence of computed concen- 
trations of intermediate-sized radicals for 
the test case for a reaction time of 0.5 ms. 
The structure for each species is presented 
in Table II. The symbols designate compu- 
tational results, the line cubic spline fits. 
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Figure 63. Temperature dependence of computed concen- 
trations of intermediate-and large-sized 
aromatic molecules for the test case at a 
reaction time of 0.5 ms. The structure 
for each species is presented in Table II. 
The symbols designate computational results, 
the lines cubic spline fits. 
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somewhere in the growth process mechanisms .h as physical 

condensation, particle coagulation and surface reactions become more 

15 *>5 

efficient and take over the soot formation process *” . In view of 
this, the model of an infinitely long reaction sequence provides only 
a lower-limit estimate of soot yield. 

As discussed in the previous section, the fused polycyclic aromatic 
molecules are formed by practically irreversible reactions. Therefore, 
Route 1 can be viewed as a sequence of similar strings of reversible 
reactions, from acenaphthylene to pyrene, from pyrene to coronene, from 
coronene to ovalene, etc,, the strings being connected by irreversible 
reactions (Fig. 64). Given such a model, the total number of carbon 
atoms accumulated in soot species can be computed exactly. In this 
computation we assumed the equilibrium constants of the corresponding 
reversible reactions in the replicating strings to be unchanged after 
coronene . 

37 38 

The mathematical formalism of this process can be presented as ’ 
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where : 


is the symbolic representation of species j within string i 
(Fig. 64); 

Pj j is coronene 

Tj ^ is the rate of formation of coronene 

X l,2 = X 3,4 = x 7,8 ~ X 10,5 = X 1 1 ,6 = k 23 lH ^ + k -24^ C 2 H ^ + k -21 
X 2,l ~ X 4,3 = x 8 , 7 = x 5 , 10 = x 6, 11 = k -23^ H 2^ + k 24^ C 2 H 2^ + k 21^ Hl 
X 2,3 = X 4,9 = X 8 , 12 = k -25^ C 2 H 2^ + k 7^ C 2 H ^ 


X 3,2 X 9,4 ~ X 1 2 , 8 k 25 fHl + k -7 


X 4,5 k 27^ C 2 H 2^ 


X 5,4 = k -27 


X 5,6 “ k 30 


X 6,5 " k -30 


X 5,9 ' k -5 


4 9 ,5 “ k 5 [H ^ 


X 6 , 7 y 8,I " k -25 fC 2 H 2 J 


and and k_£ are the rate coefficients for the forward and reverse 
directions, respectively, assigned to reaction class £ (Table IV). The 
differential eqiu Lions for the sum concentrations in this mechanism are 


* d[P ] 

2 ( LlL_ i - r 

. V dt ' r l.l 
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+ ^ X 3,2 lP i,3 ] 


® d[P. 8 1 » 


,^ { ' = ± l x X 7,8 lP i,7 1 ‘ + X 8 ,12 + 


+ 1 X 12,8^ P i,l2^ 

i=l ’ 


“ d[P ] “> 

E — h^-) = 


£ {- 
i=l 


dt 


1 x 8 f 12 lP i,8^ 


i=l 


Z x 12,8 lP i,12 ] 
i=l 


or 


dS l 

dt " r l,l " X 1,2 S 1 + X 2,1 S 2 + y 8,l S 8 


dS 2 

dt" ~ X l,2 b l " (x 2,l + X 2,3 )S 2 + X 3,2 S 3 


dS { 

dt 


=X 7,8 S 7 " tx 8,7 + X 8,12 + y 8,l )S 8 + X 12,8 S 12 


dS 


12 


dt X 8,12 S 8 " X 12,8 S 12 


where S = 2 [ p ..] , j = 1, 2,..., n 

J i=l J 


These 12 differential equations were added to the set of differential 
equations that describes the rates of change in the concentrations of 
the species which were considered individually. All mass balances were 
.aaintained exactly. Integration of the entire set of differential 
equations for the initial conditions of the test case (including Sj[t=0] 
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= 0, j = 1,2,..., 12) generated the profiles of the species concentra- 
tions and the concentration sums. 

The rat* of increase of carbon atoms in species j is 


dM. ® 


dt 


- 2 AmT . r . . + 2 Am* , . x . , . S . , 


i=l 


ij iJ 


j'*j 


J J J J J 


2 Am. . ,x. . ,S. 
j’*j JJ J 


+ 2 Am 3 :, .y., .S. , 

J J J J J 


2 Am. . , y . . , S . 

■* * ' •* ■* -i 1 -» 




JJ JJ J 


from which the total amount of soot grows as 
n dM . n oo 


dM 

dt 


= Z *5^ = j!, + j.^j C S'j K j'j + 


j = l 


2 2 (Am ,x + Am' y )S 

JJ JJ JJ J 


(5) 


where: Am. . is the measure of direct mass production from small 

reactants ; 

Am X 

jj ' is he measure of mass production associated with i,.ter- 
conversion of species with the same degree of poly- 
merization ; 

y 

AnTj , is the measure of mass production involving an increase 
in polymerization. 

In our case, one has Am X = 0 for H-atom abstractions and Aro X =Am^=2 for 
C 2 additions. The value of Am^ ^ is equal to the number of carbon atoms 
in P. . , which is 24. 

* I * 


and 


Using the above values of Am, eq.(5) takes the form 
dM 

dt = 24r l,l + 2 x 2,3 S 2 ‘ 2x ?,,2 S 3 + t2x 4,5 + 2x 4,9 )S 4 
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2x 5,4 S 5 + 2 x 6,7 S 6 + (2x 8,12 + 2y 8,l )S 8 


2x 9,4 S 9 “ 2x 12,8 d 12 


( 6 ) 


where M is the total number of carbon atoms accumulated in species 

beginning with coronene. Equation (6) with the initial condition 

M[t=0]=0 was added to the previous set of differential equations of the 

species concentrations and the concentration sums . The soot yield is 

24 

defined as in our experimental work as the number of carbon atoms 
accumulated in soot divided by the initial number of carbon atoms of 
acetylene; that is, the soot yield equals M divided by twice the initial 
concentration of acetylene. The soot growth curve using this polymeri- 
zation model is given in Fig. 65. 

As seen in Fig. 65, the model correctly predicts the existence of 
an induction time for soot appearance as observed experimentally. 
Closer examination of this figure shows that the rate of soot formation 
(the slope of the soot-yield curve) decreases after a reaction time of 
about 1 ms. This deceleration is the result of the "thermodynamic 
resistance" of the kinetic system. The reaction network is basically of 
a consecutive type and is composed of reversible reactions. Initially, 
when the "bottleneck" reactions are far from partial equilibrium, the 
dynamic development of the system exhibits the exponential growth 
typical of ''sin reactions. This is manifested in the existence of the 
induction period. As the concentrations of intermediate species grow, 
the rates of reverse reactions increase, bringing the reactions towards 
a state of partial equilibrium. When this state is attained by the 
"bottleneck" reactions, the mass flux through them decreases, which 
accounts for the deceleration in soot formation. This is probably the 
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Figure 65. A profile of soot growth for a 1.09% C 9 H--Ar 

5 3 ^ ^ 

mixture at = 3.8 x 10 mol/cm and a 
temperature of 1700 K. 
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point at which coagulation and surface growth would begin to dominate 
the process. 

Soot yields computed in the above manner at various temperatures 
are given in Fig. 66. Figure 67 presents the corresponding experimental 
results^. Except for the shift in the temperature range (see below), 
the computed soot yields compare well with the experimental ones (cf. 
Figs. 66 and 67). The model predicts both the soot bell and its shift to 
lower temperatures with reaction time. The width of the bell is in good 
agreement with experiment. The absolute values of the soot yields also 
agree well if one takes into account that the laser-extinction technique 
overestimates actual soot yields by a factor of 2 to 3^'^. Additional 
computer experiments showed that computed soot yields are dramatically 
reduced by decreases in the initirl concentration of acetylene or by 
addition of hydrogen, in accord with experimental observations. 

As for the nature of the soot-yield bell, the computational results 

24 

support our previous conjecture that competition between kinetic and 

thermodynamic factors is responsible for the phenomenon. At the lower 

temperature bound of the bell., pyrolysis reactions are rate controlling. 

As temperature increases, the thermodynamic stability of the 
41 

intermediates , particularly of small aromatic radicals, decreases and 
eventually limits the soot yield. At high temperatures fragmentation is 
dominant even at the first-ring level. The shift of the soot-yield 
maximum to lower temperatures at longer reaction times is due to the 
fact that the inflection point of soot-yield time dependence (see Fig, 
65), which is determined by the attainment of the state of partial 
equilibrium by the "bottleneck" reactions, is reached sooner at higher 



temperatures . 
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Figure 66. Computed soot yields in a 1.09% C.H.-Ar mixture 

-5 3 i L 

at p^. = 3.8 x 10 mol/cm . The symbols desig- 
nate computational results , the lines cubic 
spline fits. 
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Figure 67: Experimental soot yields in a 1.09% C-H^-Ar 

..ixture at p,. = (3.82-3.90) x 10 mol/cm 

using laser-extinction at 3.39 pm (series B of 

Ref. 24). The soot yields are calculated 

assuming Rayleigh scattering^ and using the 

39 

value of Lee and Tien for the complex 
refractive index at 3.39 pm. The symbols 
designate experimental points, the lines 
cubic spline fits. 


158 


Table VI presents the results of a sensitivity analysis for 
the reaction classes of Table IV. Somewhat surprisingly, only a small 
number of reaction classes, and sometimes only a single reaction within 
a class, have large sensitivities (Table VI). The strongest kinetic 

V 

effect, i.e, the largest value of S^, is caused by the practically 
irreversible formation of the fused polycyclic aromatics (Class 27) and 
their reactivation by hydrogen atoms (Class 23). It is interesting to 
note that the ratio £H J '/ [H^ ] is much larger than the corresponding 
equilibrium value (Figs. 54 and 59). The concentration of hydrogen 
atoms is controlled, for given initial conditions, primarily by chain 
reactions 6 and 7 (Table III) and its numerical value is well 
approximated by assuming a steady-state for these two reactions. Only 
at high temperatures, when most initial pyrolysis reactions reach a 
state of partial equilibrium, does the concentration of hydrogen atoms 
relax to its equilibrium value on the time scale of interest here. 
These results lead us to conclude that the overshoot by hydrogen atom of 
its equilibrium concentration is a driving kinetic force for soot 
formation. 

Analysis of the sensitivities with respect to equilibrium 

y 

constants, S^, indicates that the thermodynamic bottleneck for the 

growth in soot yield is at the formation of the first aromatic ring. 

Y 

The largest values of S R , obtained for Classes 13 and 8 (Table VI), are 

basically the sensitivities for the reactions 

C 2 H 2 + HC=C-CH=CH j HC=C-CH=CH-CH=CH 

and C„H_ + C.H ± HCsC-CH=CH, 

Z Z 2 *■ 


respectively. 


159 


TABLE VI. SENSITIVITIES OF SOOT YIELD WITH RESPECT TO RATE COEFFICIENTS 
AND EQUILIBRIUM CONSTANTS. 


Reaction Class 

S k 

S Y 

1 

o 
»»— • 

0 

1 

0.10 

3 

- 0.04 

- 0.20 

5 

- 0.01 

0.03 

8 

- 0.02 

0.78 

13 

0. 16 

0.79 

23 

2.02 

0.12 

24 

- 0.06 

0.06 

25 

0.27 

- 0.28 

27 

1.30 

0.24 

29 

- 0.03 

0.00 

30 

0.70 

0.56 

The sensitivities are defined^ 

as 


S l = Jen ( Y k / Y 0 )/£ n ( l / 5 ) 

and 

= i!n(Y K /Y 0 )/£n5, 


where Y Q is the soot yield computed for 1.09% C 2 H 2 -Ar mixture at = 
3.80 x 10 5 mol/cm 3 , T,. = 1700K, and reaction time of 0.5 ras ; Y^ and Y^ 
are the soot yields computed for the same conditions but with the rate 
coefficient reduced by a factor of 5 and with the equilibrium constant 
increased by a factor of 5, respectively, for a given reaction class. 

The reaction classes omitted have absolute sensitivity values less 




than 0.005. 
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The results of the sensitivity analysis (Table VI) indicate a 
strong interaction between rate coefficients and equilibrium constants. 
This we believe to be the cause for the mismatch in the temperature 
ranges of the experimental and numerically simulated soot-yield bells. 
Further study of the thermochemical and rate constant data and the 
mechanism itself should reveal the cause of this mismatch. 
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MAIN CONCLUSIONS 

1. The presented experimental results futher support our earlier 
conjecture that the soot formation mechanism is the same for both 
pyrolysis and oxidation of hydrocarbons. 

2. The experimental results performed to investigate the effect of 

fuel structure on soot formation can be rationalized within the 
framework of our conceptual models and thus support them. One 
deviation is the results obtained with vinylacetylene . These 

results indicate that conjugation by itself is not a sufficient 
factor for determining the sooting tendency of a molecule - 
structural reactivity in the context of the chemical kinetics is 
the dominant factor in soot formation. The latter conclusion is 
also supported by the results of our detailed kinetic modeling. 

3. It was demonstrated that soot formation in shock-tube pyrolysis of 

acetylene can be explained by a detailed ki. etic mechanism composed 
of conventional hydrocarbon reactions. The main mass growth was 
found to proceed through a single dominant route in which fused 
polycyclic aromatics play a particularly important role: tdeir 

practically irreversible formation reactions have the effect of 
"pulling" chains of reversible reactions in which H-atoms 

reactivate aromatic molecules. The main reaction bottleneck appears 
at the formation of the first aromatic ring. The overshoot by 
hydrogen atom of its equilibrium concentration is a driving 

kinetic force for soot formation. The model is in accord with 

product distributions observed in flames. 


fZ 
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FUTURE RESEARCH DIRECTIONS 

The ultimate goal of this program is the development of a detailed 
kinetic model for soot formation in combustion of hydrocarbon fuels. To 
achieve this goal, and taking into account our present results, we will 
concentrate on the following objectives: 

1. The continuation of the experimental investigation of soot 
formation in fuel mixtures. 

2. Refinement of the developed chemical kinetic model of soot 
formation in acetylene pyrolysis. 

3. Extension of the developed chemical kinetic model of soot formation 
to other conditions. Of immediate interest are acetylene oxidation 
and pyrolysis of other representative fuels. 
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